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Abstract 

Data obtained from a spatial and temporal distri- 
bution study of sewage effluent, traced with fluorescent 
dye, discharged from the Municipality of Metropolitan 
Seattle (METRO) West Point treatment facility on Puget 
Sound are presented. 

Following an initial dilution by a factor of about 
l/lOO at the discharge line diffuser ports, effluent 
diluted further by only 1/10 was found approximately 
6000 meters north (ebbing current) and 3500 meters south 
(flooding current). Significant concentrations diluted 
by a factor of about 1/3 after leaving the diffuser ports 
were found near shore, and concentrations diluted by a 
factor of 1/2.3 were found on the West Point beach (3.7 
ppt effluent). 

From contours of sewage concentration and density, 
wide spatial and temporal variations were observed in 
both throughout the 5 days that the experiments were con- 
ducted '30 August through 3 September 1974-). Transverse 
cuts of the plumes indicated that contours of sewage ef- 
fluent generally coincide with density contours indi- 
cating that density stratification is influential in 
confining the distribution of sewage in the Sound. A 
longitudinal cut, however, indicated significant pene- 
tration of density layers by the effluent near the ends 
of the effluent plume. Implications of these phenomenon 
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on modeling effluent diffusion in the Sound are discussed. 
A general "buildup of dye in the Sound over the period of 
the experiments was not observed. 

Comparison of these dye studies with drogue studies 
carried out simultaneously show good agreement quanti- 
tatively for the 30 August data and qualitatively for 
data obtained on other days. 



I. Introduction 



From 30 August to 3 September, 1974, the Applied 
Physics Laboratory (APL) of the University of Washington 
undertook a study of the sewage effluent distribution in 
Puget. Sound discharged by the Municipality of Metropoli- 
tan Seattle (METRO) processing plant located on West Point, 
Seattle. This study addresses two principal questions: 
Where does the sewage go? What is the nature of the mix- 
ing process acting on the sewage as it is transported 
through^ Puget Sound? By observing the nature of the ans- 
wers to these two questions, conclusions regarding the 
applicability of model types can be drawn. Whether a 
buildup of sewage concentration occurs or whether sewage 
residence times at detectable concentrations are limited 
is investigated as well. 

Data were obtained by cycling a towed sensor, "fish," 
through depths where the sewage was detected. Fluores- 
cent dye was added to the sewage effluent as it was dis- 
charged from the processing plant. This dye was detected 
by the fish *as it was towed through the effluent in the 
Sound. On board, real-time instrumentation allowed ac- 
curate determination of dye concentration, fish depth, 
and vessel position with respect to the sewage effluent 
discharge ports (outfall). The effluent was traced from 
the outfall to the areas in the Sound that it reached. 

The fish was capable of determining temperature and 
conductivity. Density data determined from the tempera- 
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ture and conductivity data are also used in this analysis. 

Dye concentration, density, depth and vessel loca- 
tion data were recorded on magnetic tape on board and then 
processed ashore to provide raw information used to hand 
contour the data. The visualizations permitted by the 
contoured data lead to most of the observations made in 
this report. 

An investigation of current properties using drogue 
information was conducted by Evans -Hamilton, Inc. in paral- 
lel with the dye studies.'®' Comments concerning the con- 
clusions reached by these studies are made as a result of 
the data obtained from the dye studies. 

This investigator prepared several programs for use 
with the hand held computer-calculator, HP-65, for navi- 
gation and data work-up that are presented in Appendix H. 

p 

II . Experimental Procedure 

Rhodamine B dye was injected into the sewage efflu- 
ent just upstream of the post-chlorinators at a constant 
rate of 1.73 kg/hr (dry weight) fo^ the first four days 
of the experiment (30 August - 2 September). The flow 
from the outfall (located about 1100 meters west of West 
Point at a depth of 73 meters) was 70.9 million gallons/ 
day (mgd) over this period yielding an average concentra- 
tion of dye in the pipe of 1.5 x 10” g/cc. On 30 August 
the average flow was 82.8 mgd and the dye injection rate 
was 1.85 kg/hr yielding a concentration in the outfall 
of 1.4 x 10”^ g/cc. A dye injection rate six times that 
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for the first four days was planned for the fifth day, 
but cancelled due to towing cable failure. 

Dye concentration was measured with an instrument 
package ( "fish" ) towed behind APL 1 s 50-foot research boat 
and cycled over a depth appropriate for detecting the dye. 
Prior to pumping dye into the system, a background fluor- 
escence level equivalent to 5 x lO - ^ g/cc dye concentra- 
tion was determined in Puget Sound waters. Since dye 
concentration in the outfall is about 1.5 x 10~^ g/cc, 
this means that a minimum effluent concentration of 0.33 
parts per thousand (ppt) could be observed. This cor- 
responds to a maximum dilution factor of 1/3000 obser- 
vable. Towing speeds varied from 3 to 6 knots. The fish 
consists of temperature, electrical conductivity, dye con- 
centration, and depth (pressure) sensors housed in a 
streamlined aluminum tube. A pump within the fish forces 
water through the sensors. An electronics package accepts 
sensor output, multiplexes the data onto a single conduc- 
tor within the towing cable and transmits the information 
to shipboard equipment. Figure 1 shows the general con- 
figuration of the fish and internal equipment. A hydrau- 
lic winch was used to tow and cycle the fish vertically. 
This gasoline powered unit is capable of hauling in at a 
rate of 1 meter/second. Both water depth and fish depth 
indications were displayed at the operators station 
enabling the operator to cycle the fish near the bottom 
while avoiding a collision with it. 
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The boat's position was determined to within a few 
meters using a Motorola Mini-Ranger radar navigation 
system - ^, whi ch determines range between each of two on- 
shore transponder units and the shipboard master trans- 
mitting unit. Ranges were continuously displayed on the 
face of the master console aboard the boat. An HP-65 pro- 
grammable calculator was used to convert range-range in- 
formation to X-Y coordinates based on the sewage outfall 
position. This real time navigational information was 
essential to assure that the boat's navigation closely 
followed desired paths with respect to the outfall. 

Ranges, dye concentration, fish and water depths, 
and clock time were ail available aboard during the con- 
duct of the experiment. Range, all fish data, water depth, 
time, and manually entered "event number" were recorded 
once per second on magnetic tape that is compatible with 
APL's CHI 2130 computer. 

Prior to the experiment, a series of observations 
was made of the UW Oceanography Department's Puget Sound 
Model in order to provide a rough feeling for the expec- 
ted behavior of the dye. A wide degree of variability from 
tidal cycle to tidal cycle suggested that a preset pattern 
for sampling would not be workable. Instead, it was 
decided that a series of arcs crossing the Sound in east- 
west directions starting near the outfall and working 
away in the direction of the current would be the most 
productive means for detecting the area containing the 
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effluent. In addition to follow ins this general stra- 
tegy, sweeps hack up the plume (s) were conducted; and 
several excursions through water outside the boundaries 
of well defined plumes were made to detect detached 
pools that observations of the Puget Sound Model suggested 
might exist. 

These studies and the drogue studies were coordina- 
ted in that the depth at which dye was detected was com- 
municated to Evans-Hamilton personnel for use in setting 
drogue depth, and the establishment of flow following 
slack water was determined from the drogues and used to 
time the commencement of the dye detection sweeps, nor- 
mally, about 30 minutes were allowed to pa.ss after flow 
commenced in order to assure a well established current 
pattern prior to sampling. 

III. Dye Concentration and Density Contouring 

Data obtained aboard on magnetic tapes was returned 
to the APL and processed on the CHI 2130 computer to plot 
dye concentration and density profiles for each cycle of 
the fish. Figure 2 is a sample of this output that was 
used in plotting dye and density contours. The boat's 
position at the beginning and end of each profile was 
determined from a printout of the data in tabular form 
as shown in Figure 3. Using this information, plots can 
be made that correspond to the X and Y positions the 
ship traversed. The tracks made during the dye detection 
sweeps, however, closely approximated straight lines. 
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FIGURE 3 - TABULAR DATA OUTPUT 
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This allowed simplifying the contoured data to he plotted 
with respect to distance along the approximate straight 
line track. The straight line tracks are shown as the 
first figure of each of Appendices A through F. Each 
plot that follows in these Appendices shows distance 
along the respective tracks (corresponding to the ap- 
propriate event numbers) with the z ero point correspond- 
ing to X = 0 or Y = 0 depending upon whether the slope 
of the track lies within 45° of the X or the-Y axis res- 
pectively. Where density contours have heen plotted, 
they were done on translucent paper so that they may he 
overlaid on the corresponding dye -contours . This method 
greatly improved the visualization of the contoured data 
and leads to some of the conclusions expressed in this 
report. The depth dimension used in each case is 10 
meters/inch with the water surface at the top line. All 
contoured data are presented as though the observer were 
looking north for the east-west cuts across the plumes, 
and as though the observer were looking west for the north 
south cuts. * Each data point on both the dye concentra- 
tion and density contours is indicated by a symbol that 
represents the observed concentration at the point plotted 
Figure 4 provides a listing of all symbols used. 

The reader will note upon a cursory inspection of 
the plotted data that it is closely spaced along the 
path of the fi3h compared to the horizontal distance bet- 
ween successive depth .cycles. During the preliminary 
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Figure 4 

Contour Plot Symbols 
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stages of preparing the data, computer programs that were 
capable of contouring were investigated and it was found 
that none were capable of handling data that was so 
unevenly distributed. Therefore, hand contouring became 
necessary. Although, in most cases, contours could be 
drawn in only one way, occasionally, the data became 
ambiguous and a judgement regarding the most appropriate 
way to draw the contour had to be made. 

The contoured data are presented in six. segments 
in Appendices A through F. These correspond to the data 
obtained on the sweeps indicated in Table 1. 

(ft 

Table 1 

Experiment Descriptions 



Day 


Date 


Flood/ 

Ebb 


Forth/ 
South of 
Outfall 


Direction of 
Sweep (Toward 
or Away from 
Outfall) 


Appendix 


1 


30 Aug. 


Flood 


South 


Away 


A 


2 


31 Aug. 


Ebb 


North 


Away* 


3 


2 


31 Aug. 


Flood 


South 


Away 


C 


2 


31 Aug; 


Flood 


South 


Toward 


D 


4 


2 Sep. 


Ebb 


North 


Single loop 
made in area 


2 


5 


3 Sep. 


Ebb 


North 


Single loop 
made in area 


F 



The last event of this segment was a sweep back up the 
plume toward the outfall. 



The time each new event was begun is indicated on 
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each track segment on the track plots. The reader will 
note the significant time spans encompassed in the pat- 
terns presented. It is obvious that the data presented 
is not synoptic. Over the period encompassed during 
single events, the contours probably do approximate a 
synoptic - picture in that during the short period of time 
(about 30 minutes more or less) the pattern existing in 
the Sound is unlikely to be changed very much. Further- 
more, changes that do take place within this time frame 
are likely to conform to the random processes which pro- 
duce a given pattern that exists at any given instant. 

The sane statement is not observed to apply over time 
scales commensurate with that of the tidal cycle itself. 
Examples of significant change in the plume patterns 
are found in the data during lamer sweeps of the plume 
area within the same tidal C3 r cle. 

The times of the slack waters that preceded and fol- 
lowed the plotted tracks, as determined from tidal current 
tables, are also noted on each track plot. 

IV. Observations 

The data, as presented in Appendices A through F, 
are intended to be useful to the reader for the purposes 
of whatever perspective might be appropriate to his in- 
terests. Pertaining to the specific objectives of the 
METRO study, however, some observations are offered. 
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Variability of Sewage Concentration and Density 

The immediate and obvious observation upon inspec- 
tion of the contoured data is the spatial variability of 
both dye concentration and density that was encountered. 

Of all the data obtained, day 1 contours are the most 
uniform,' but even these show significant variations in 
depth and along the track. Some of these events show; 
two areas of high dye concentration instead of a single, 
uniformly developed peak (see events 7 and 10.2, day 1) 
that might be expected of an "idealized" plume. Intru- 
sions of clean water (concentration of dye less than 
background) are seen to divide the plume in depth also 
(see events 11, 12, and 13, day 1). 

The data obtained on days 2, 4, and 5 all show pat- 
terns of dye concentration more complex than that of day 
1. Four or more discrete filaments are discernible in 
much of this data, providing an indication of the complex 
mixing processes occuring in the Sound. It bears mention- 
ing, again, that caution must be exercised in analyzing 
the contoured data presented since a significant amount 
of time passes from the start to the completion of each 
event. 

Models and Implications 

The foregoing observations bear significant impact 
on the types of models that may be devised to describe 
the effluent distribution in the Sound. Many models 
assume that the diffusion is radially symmetric, and, 
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further, that probability distributions are Gaussian or 
of a shape similar to that of a Gaussian distribution. 

From the foregoing observations, it is evident that nei- 
ther of these conditions hold in the Sound. The day 1 
data most closely approximates these conditions, but 
even this data implies that such models must be used in 
a very limited fashion. It has been noted, however , that 
although the distributions observed hardly conform to 
the conditions mentioned above, a good fit was obtained 
on day 1 for a model that assumes the "centerline" con- 
centration, to be approximately an inverse function of 
distance from the outfal 1 (concentration proportional to 
l/Y a , where 1 <C a 3/2 ) 4 . The data from day 4 appear to 
conform to such a model as well, but only 3 events were 
obtained across the plume(s), and the data is, therefore, 
not conclusive. Certainly, any model that attempts to 
span a time scale greater than one tide could not suc- 
cessfully incorporate the assumptions of a smooth dis- 
tribution and radial symmetry. 

A further assumption, that the effluent may be des- 
cribed in two horizontal dimensions, is frequently made 
in models used to describe flows such as under study here. 

In order for such an assumption to be valid, sufficient 
stratification must be present to confine the effluent 
flow along density contours only. It will be seen in the 
following observations that, although stratification greatly 
affects the distribution of effluent, this assumption, 
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also, is not valid during portions of the tidal cycle. 

Depth Dis tr ibution of Ef fl uent 

From the overlays of density contours it is seen 
that the stratification in the Sound is influential in 
the observed distribution of sewage as indicated by the 
dye. Little evidence of dye contours crossing density 
contours exists in the east-west crossings of the plume(s) 
on day 1. Subsequently, on day 2 and 4 data for instance, 
the trend for dye contours to follow density remains evi- 
dent, although frequent crossings do occur. This may 
be partially due to the method in which the data was 
processed. Following the first day's experiment, a sig- 
nificant lag in density detec tor (s) response occurred, 
possibly due to a pump or piping malfunction (leak) that 
went undetected. This necessitated averaging the density 
data between successive half cycles to eliminate the time 
lag. Therefore, the contours presented for days 2 and 4 
represent values averaged over from 50 to 200 meters, 
depending upon the distance covered by the boat between 
cycles. .The smaller scale density fluctuations were 
thus averaged out while the dye detection equipment con- 
tinued to function properly, allowing the smaller scale 
dye concentration structure to be observed. 

It is noteworthy that, in areas very close to the 
outfall and at low current velocities, the effluent was 
able to significantly affect the density structure. 

Event 5 on day 1, event 12.1 on day 2, and event 2 on 
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day 4 all show a divergence of density contours surround- 
ing areas of high dye (sewage) concentration. These areas 
may be considered "active" in the sense that the deeper 
portions of the high sewage concentrations are lighter 
than their surroundings and the shallower portions are 
heavier.- These conditions were generally observed within 
about 400 meters or less of the outfall. During periods 
near the outfall when the current velocity is relatively 
low, significant vertical penetration of the dye appears 
to take place from these active processes and may lead to 
the observations described in the following paragraphs. 

Longitudinal Distribution o f Sewage 

Event 11 on day 2 is a longitudinal track from north 



to south near the end of an ebbing current. The density 
contours overlaid upon the dye contours clearly show that 
dye penetrates density contours. The plume, shown to be 
separated by about 1400 meters in the portion generated 
during the latter half of the tidal cycle, displays ends 
in each segment that reach a lesser density layer at the 
northern end than is the case with the remainder of the 
segment. Because the density contours are penetrated, 
models that attempt a two dimensional description cannot 
be directly applied. (It has already been shown that a 
simple plume description is not appropriate due to the 
extreme variability of the flow both spatially and tem- 
porally. ) The cause for this observed penetration of 
the density layers cannot be stated with much confi- 
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dence "because of the limited amount of data available 
longitudinally. Because shallow portions of the plume 
segments lie at the ends, it is possible that the greater 
penetration occurs during periods of low current vel- 
ocity that occurs during a change of tide and due to 
changing current conditions in the middle of a tidal 
cycle (perhaps the disestablishment of an inshore eddy 
significantly slows the current flow over the outfall 
for a period of time followed by a return to. greater 
velocities as the new trajectory becomes established). 
During this period, lateral mixing processes might be more 
severely restricted and thereby diffuse less of the 



buoyant momentum of the effluent as it leaves th 
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ports of the outfall. This would tend to cause the ef- 
fluent to penetrate to a shallower density area than if 
the current velocity were higher. Further investigation 
is needed to better clarify the processes that cause this 
phenomenon. 

Surface Penetration 

It is significant to note thau no dye was detected 
on the surface during these experiments. In general, 
only infrequently did dye penetrate to above 20 meters, 
reaching only about 4 or 5 meters at the most. A signi- 
ficant exception to this occurred on day 5, from 2000 
to 3500 meters north of the outfall. Significant con- 
centrations were detected to about 4 meters of the surface, 
and it is likely from .observations of the contoured data 
of events 0 and 2.2 of day 5, that dye may have reached 
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a depth of 1 or 2 meters. Significantly, wind condi- 
tions that prevailed during day 5 were different (though 
not particularly severe) than on the previous days that 
the experiments were conducted. Wind records from the 
West Point Coast Guard station show a reversal of the wind 
from NS on day 4 to SSS on day 5 with an increase in 
average speed from 7.5 to 10.7 knots. Although the mag- 
nitude of this speed increase does not seem large, a 
significant change of surface wave conditions were noted 
on day 5. Also, a lesser density stratification was noted. 
These conditions obviously had an impact on the dye dis- 
tribution that was observed. It would be expected that 
an even greater distribution in depth might be noted In 
the winter, and this, from preliminary indications thus 
far in the winter experiment analysis, was the case. 

Dilution of Sewage Effluent in the Sound 

The maximum concentration of dye observed through- 
out the conduct of the experiments occurred during a 
flood tide on day 2, event 12.1. A value of 1.26 x 10~^ g/cc 
was observed 400 meters south and 200 meters east of the 
outfall. This represents a dilution factor of 1/119 
since the mean concentration in the outfall on day 2 
was 1.5 x 10“ g/cc. Of interest is that a concentra- 
tion of 1.6 x 10“ 10 g/cc was observed over 6000 meters 
north of the outfall on event 10, day 2. This represents 
a dilution factor of only 1/7.9 in the distance.' Since 
the subsequent event, 11 on day 2, shows two segments, 
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the impression that the northerly segment belongs to 
a previous tidal cycle was gained, but found not likely. 
Direct measurement of current velocity by drogues was 
not accomplished on this day's experiment. However, 
by using the maximum current observed for all of the 
drogue studies (day 1), of 38 cm/sec, and the fact that 
the average current during a tidal current of the dur- 
ation observed is approximately .64 times the maximum 
velocity, a current velocity average much in excess of 
24.2 cm/sec over the tidal cycle involved is not likely. 

In six hours at this speed, a distance of 3230 meters 
could have been traveled by the leading edge of the 
plume without taking into consideration mixing that 
occurs. It is not unreasonable, therefore, that it 
reached 6000 meters in a single tidal ebb current. Pos- 
sibly, during the portion of the ebb that occurred at 
relatively high current velocities, the plume was directed 
toward the eastern part of the Sound (toward Shilshole Bay) 
and as the ebb began and also near (but not at) the 
end when it slowed, it changed to a more central route 
up the center of the Sound. Such a pattern (and the many 
others that represent such widely variable character- 
istics) further serves to emphasize the difficulty in 
modeling flow in the Sound. 

From these observations, it is evident that, as 
pointed out in the preliminary report to METRO, most of 
the mixing and dilution of the effluent occurs at the 
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diffuser ports at the outfall. High concentrations 
we re observed south of the outfall as well, although at 
a distance only half that observed to the north. On day 
2, event 20, a concentration of 4.0 x lO - ^ g/cc was 
observed nearly 3000 meters to the south and very near 
shore. This concentration represents a dilution factor 
of only 1/3.2 to the highest concentration observed near 
the outfall on this day. 

The data indicates that sewage reaches considerable 
distances north and south of the outfall and very near 
shore with little dilution from mixing in the Sound, 
within one tidal cycle. 

Buildup of Sewage 

Throughout the five days the experiments v, r ere con- 
ducted, a buildup of a dye background was not observed. 
All dye observed showed distinct boundaries, whether 
part of the main plune(s) or drifting patches. The con- 
centrations within the plume areas remained relatively 
consistent and clean (less than background dye concen- 
tration) water can be found bounding and intruding within 
the main concentrations. Sever 3.1 excursions were made 
outside the main plume(s) and dye was not detected. It 
is difficult to positively identify patches of dye as 
belonging to a previous tidal cycle, although event 4 
on day 2 is a likely candidate since the location and 
concentration of this dye patch seems to indicate that 
it is not part of the main plume(s). It appears that 
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most of the sewage that enters the Sound may he aiixed 
to below background concentrations within two tidal cycles. 

Sewage Hear Shore 

Concentrations of effluent ranging as high as 2.7 
ppt were observed on several occasions very near shore, 
both north and south of the outfall. This is somewhat 
evident during the day 1 experiment in events 7 and 11, 
but most vivid on day 2 during events 20 and 21 taken 
during southerly sweeps. Observations during northerly 
sweeps include day 4 event 3.1> and day 5, events 0 and 2.2, 
As shore is neared, water depth decreases and the depth 
to which the fish is cycled is correspondingly reduced 
and this is seen in the plots of the above tracks. Where 
the fish depth decreases in these plots significant con- 
centrations of dye were detected, and, in most cases, 
the concentrations actually increase, possibly due to a 
near shore eddy. That effluent reaches shore is confirmed 
by bottle samples taken on both the north and south beach- 
es of West Point. Dye concentrations as high as 3.1 
x 10~ 1G g/cc (corresponding to an effluent concentration 
of 3.4 ppt) were observed on the beach north of West 
Point and as high as 5.5 x 10 - ^" G g/cc (3.7 ppt effluent) 
south. It is evident that sewage diluted by a factor of 

1/2.3 of that found to be the maximum after exiting the 

♦ 

diffuser ports finds its way ashore. 



V. Amount of D' r e in the System 



In an attempt to determine whether significant por- 
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tions of the effluent discharged into the Sound were 
going undetected, a mass balance of day 1 data for events 
10.2 through 12 was made. Day 1 was chosen since the plume 
pattern was more uniform that day and a good cross sec- 
tion of the dye in the plume was obtained for each of the 
events used. Dye contoured data was averaged for each 
5 m by 25 m area transverse to the plume. These values 
were then added to provide a figure for the total amount 
of dye per longitudinal centimeter as demonstrated by 
the following relation: 



C ev [Vim]= (£ CtLl^Mo+cwbl) 25 x io 4 [cm 2 J n 

L" I 



where: Cj_ = the concentration observed in an indi- 



vidual 1.25 x 10^ cm^ cross sectional 
area (g/cm^/l.25 x 10^ cm^) 
c ev = a linear concentration along the plume 
axis (g/cm) for a given event (ev) 

6 p 

n = the number of 1.25 x 10 cm squares 
for which a concentration (C^) was 
observed during the event 



The values of c ev for the three 'events were averaged, 
thereby entertaining the approximation that current pro- 
perties were the same as each of these segments of water 
passed the outfall. Using the average value of c ev = 
3.08 x 10” (g/cm) , and assuming that this average fig- 
ure persisted for one. hour (t = 5600 sec), and further 
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assuming that this segment of dye passed over the out- 
fall at an average of one hour following the commence- 
ment of the flood tide (an approximation based on the time 
of placement of the drogues into the water and the observed 
position of the boat with respect to the drogues during 
this portion of the experiment) at a velocity half L of 
the maximum observed from the drogue data (v = .5 x 38 
cm/sec), the quantity of dye, q g/hr, entering the Sound 
is determined as follows: 

q = c ev x v x t 

= (3.08 x 10~ 2 )(0.5 x 38 ) (3600) 

= 2100 g/hr (or 2.1 kg/hr) 

The actual injection rate measured at the dye pumping 
equipment located at the West Point sewage treatment 
facility was 1.73 kg/hr. Agreement between these two 
figures is good and there is no indication that dye is 
entering the system and going undetected into areas 
outside the observed plume. 

VI . Comparisons with Drogue Studies . 

The data obtained on day 1 were digitized and ana- 
lyzed to provide statistics that could be directly com- 
pared with those obtained by C. C. Ebbesmeyer of Evans- 
Hamilton, Inc. Density contours for events 6 through 14 
of day one were used to prepare digitized data for each 
5 m by 25 m square cross-section (1.25 x 10 cm^). 



These 
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data are presented in Appendix G. Values shewn are in 
tenths of a power of ten that is added to the power -10.3 
of ten. For example; a value of 4.5 shown on the grid is 
equivalent to a power of 0.45 to be added to -10.3, 
yielding a power of ten equal to -9.85. The concentra- 
tion represented is therefore 10”^*^^ in this case, or 

_Q g O 

0.14 x 10 over the 1.25 x 10 cm area. The abscissa 
dimensions are arbitrarily selected digits each repre- 
senting 25 meters along the track. A program was pre- 
pared for the HP-65 computer-calculator to calculate 
centroids in both dimensions, and another prepared to 
calculate the standard deviations .and variances of the 
distributions. These programs are included as Appendix H. 

The data determined in this way are presented as 
Table 2. A plot of a plume that corresponds to Figure 
3a of the Evans-Hamilton final drogue report^ is pre- 
sented as Figure 5. Standard deviation data were converted 
to represent those obtained by C. C. Ebbesmeyer along 
the axis shown in his Figure 7a. A plot of this data 
versus distance from the outfall is presented as Figure 6, 
and can be directly compared with Figure 8 of the Evans- 
Hamilton report that is enclosed' as Figure 7 (overlay 
Figure 6 on 7 for comparison). 

The day 1 dye experiment analyzed in this way gen- 



erally confirms the drogue data. Differences that do 
occur are probably more a result of the difference in 
time that the dye data was obtained with respect to that 
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Table 2 









Day 1 Dye 


Statistics 






Event 
Ho . 


X 

(m) 


V 

(m) 


(lcA cm) 


(10 2 '' cm) 


O 

X 

CO 

O 

i — 1 


(10 6 ^cm) 


6 


•188 


53.5 


1.99 


10.8 


. 2.96 


1.18 


7 


183 


57.2 


2.21 


12.7 


4.91 


1.62 


10(1) 


160 


50.6 


2.49 


10.0 


6.21 


1.00 


10(2) 


359 


54.0 


2.65 


15.0 


7.03 


2.25 


11 


188 


57.3 


2.42 


14.7 


5.87 


2.17 


12 


294 


54.9 


2.40 


16.2 


5.78 


2.61 


13 


411 


59.6 


3.78 


14.9 


14.31 


2.21 


14 


315 


56.4 


4.77 


13.9 * 


22.73 


1.93 



’ V 



^7 





Figure 8. Average .standard deviatic/n ( x4 ) of simulated plumes versus 
distance from outfallpiQURE 6 
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Figure 8. Average .standard deviation ( x4 ) of simulated plumes versus 
distance from outfall. 



Figure 7 
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of the drogue data (dye data is about to 1 hour later 
than the drogue data oh day 1). Also, a qualitative 
inspection of data from experiments conducted on subse- 
quent days confirms the general pattern observed using 
the drogues. However, from the dye studies, it is seen 
that distinct filaments are observed along with a wide 
degree of variability of concentration in general, and 
these observations would indicate that the current pro- 
perties as observed using drogue data are likely to be 
significantly optimistic with respect to the concentra- 
tions that _ they predict. It is not possible to detect 
individual filaments using drogues, and it is evident 



hat multi -filament "plumes" persist throughout the data 



th 



Eddy coefficients therefore cannot generally be predicted 
confidently using models consistent with the smoothly 
averaged type of data obtainable from drogues. 



VII . Summary 

Most of the mixing over a single tidal cycle ap- 
parently takes place at the diffuser port at the outfall. 
The maximum concentration found in the vicinity of the 
outfall was 1.26 x 10**^ g/cc, representing a dilution of 
the sewage by a factor of 1/119. After leaving the vici- 
nity of the outfall, very large spatial and temporal 
variations of sewage concentrations are found in the 
Sound. (The same can be said of density profiles.) In- 
dications of sewage diluted by a factor of less than 
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1/10 in the Sound proper were found as far as 6 km north 
and 3.5 kin south of the outfall. Near shore, concentra- 
tions as high as 4 x lO - ' 1 ' 0 g/cc were noted, representing 
a dilution by a factor of 1/3.1 since entering the Sound. 
Sample bottles ashore were found to contain concentra- 
tions as high as 5.5 x 10 -10 g/cc, representing a dilution 
by a factor of just 1/2.3. 

The data indicate that the density stratification is 
influential in confining the distribution of sewage in 
depth, and little significant crossing of dye and den- 
sity contours occurs in the transverse cuts across the 
plunes. A longitudinal cut back up a plume near the end 
of an ebb current, however, demonstrated significant cross- 
ing of density contours by dye concentration contours 
which was most pronounced near the ends of the plumes. 

This may indicate that the sewage is more effective in 
penetrating the pycnocline over the outfall during low 
current velocities. 

Models that are used to describe diffusion of the 
nature occuring from the METRO sewage outfall at best 
Point often assume that the flow is two dimensional, or 
even that the diffusion is radially symmetric. Experi- 
mental observations indicate that none of these assumptions 
can be applied to the conditions existing in the Sound 
without accepting very significant error except over an 
extremely limited range of space and time scales. If 
two-dimensional models are employed, they should be ap- 



31 



plied to that portion of the tidal cycles that contain 
well established currents, since some evidence exists 
that sewage penetrates density contours more during the 
initial stages of the cycle. A model that assumes the 
"centerline" concentration to be approximately an inverse 
function- of distance from the outfall was applied to 
the data with limited success. Such a model applies 
during tidal cycles that produce a plume that is well 
defined and contains little transverse variation. Such 
a condition leads to less mixing than in plumes contain- 
ing more variability (ie: more filaments or eddies). 

Only the data from day 1 provided a significant indica- 
tion for success of such a model. 

The buildup of a general dye background was not ob- 
served although it is not possible to say that the clear 
water observed in the experiment was a result of advec- 
tion of dye out of the observation range (unlikely), or 
a result of dye having been diluted to below the limits 
of detectability. This experiment did not address the 
question of what happens after a tidal cycle is complete 
with, respect to whether dye is transported back across 
the outfall, thereby surviving the current reversal. 

The winter experiment, under analysis at this time, in- 
dicates that significant patches of dye do, in fact, 
survive the reversal in current and are transported back 

Q 

across the outfall. 

Comparison of the dye studies data with drogue stu- 
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dies data indicates that the two are generally consistent. 
Significantly , however, the variability of conditions 
that were found in the dye studies cannot be determined 
by drogue studies. Therefore, the drogue data probably 
underestimates the peak concentrations likely to be 



observed • in nature . 



33 



References 



1. C. C. Ebbesmeyer and Akira Okubo, A Study of Current 
Properties and Mixing Using Drogue Movements Observed 
Hear West Point in Puget Sound, Washington , Ev an s - 
Hamilton, Inc., Seattle, Washington, (November 1974), 
report prepared for the Municipality of Metropolitan 
Seattle, Seattle, Washington. 

2. William ?. Bendiner, Dispersion of Effluent from the 
West Point Outfall Interim Report , Applied Physics 
Laboratory, University of Washington, (15 January 
1975), report prepared for the Municipality of Metro- 
politan Seattle, Seattle, Washington, pp. 1-4. 

The experimental procedure outlined in this re- 
port closely parallels the referenced document. 

3. Detailed information concerning the Motorola Mini- 
Ranger III system used in this experiment may be 
obtained from Mr. Ernest H. Linger of the Applied 
Physics Laboratory, University of Washington, Seattle. 

4. Bendiner, op. cit., p. 5. 

5. Ibid., p. 6. 

6. U. S. Department of Commerce, National Oceanic and 
Atmospheric Administration, Tidal Current Tables 1974 , 
Rockville, Kd.: National Ocean Survey, 1973, P. 232. 

Water passing over the outfall one hour after 
the commencement of the flood tide possessing an in- 
terval between slack and maximum current of 3 hours 
has a velocity .5 of the maximum current according 
to Table A. 

7. Ebbesmeyer, op. cit., p. 5. 

8. Personal communication with Mr. William P. Bendiner 
of the Applied Physics Laboratory, University of 
Washington, Seattle, March 31, 1975. 



34 

Appendix A 



Day 
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Flood/ 

Eb’o 


North/ 
South of 
Outfall 


Direction of 
Sweep ( Toward 
or Away from 
Outfall) 


1 


30 Aug. 


Flood 


South 


Away 



Event No. 


Type of Contour 
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(Track Chart 




, . 35) 


05 


Density 
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Dye 
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50 
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14 


Dye 
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Outfall) 
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* The last event of this segment was a' sweep hack up 
the plume toward the outfall. 
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Type of Contour 
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(Track Chart 
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Day 


Date 
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North/ 
South of 
Outfall 


Direction of 
Sweep (Toward 
or Away from 
Outfall ) 
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31 Aug. 
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South 


Toward 



Event Ho. 


Type of Contour 


Page Ho. 


(Track Chart 




, . 74) 


20 


Dye 


75 


21 


Dye 


76 


22 


Dye 
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Day 
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Flood/ 
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Direction of 
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South of 


Sweep ( Toward 








Outfall 


or Away from 
Outfall 


4 


2 Sep. 


Ebb 


North 


Single loop 
made in area 



Event No. 


Type of Contour 


Page No. 


(Track Chart 


. . 79) 


1 


Dye 


80 


2 _ 


Density 


81 


2 
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82 


3.1 


Dye 




3.2 


Dye 


84 


4.2 
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South of 
Outfall 


Direction of 
Sweep (Toward 
or Away from 
Outfall 


5 


3 Sep. 


Ebb 


North 
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made in area 



Event No. 


Type of Contour 


Page Ho. 


(Track Chart 
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2.1 
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Appendix G 

Raw Statistical Data 
Day 1 



Event 


Page 


06 


92 


07 


93 


10(1) 


94 


10(2) 


95 


11 


96 


12 


97 


. 13 


93 


14 


99 



100 



Appendix H 
HP -65 Programs 

Title Page 

METRO - Navigation Program 101 

Nav. Prograra Trigonometric 

diagram 102 

Centroids (Dye Data) 103 

Variance and Standard Devi- 
ation (Dye Data) 104 

Two Bearing Fix 105 

Sextant Conversion to True 

Bearing 106 
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